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We present an SU (4) model of high-T c superconductivity. One dynamical symmetry of this model 
corresponds to the previously proposed SO(5) model for unification of superconductivity and an- 
tiferromagnetism, but there are two additional dynamical symmetries: S l O(4), associated with an- 
tiferromagnetic order and SU(2), associated with a D-wave pairing condensate. These provide a 
3-phase microscopic model of high-T c superconductivity and permit a clear understanding of the 
role played by the 5*0(5) symmetry. 



There are strong arguments that the mechanism for 
high-temperature superconductivity (SC) is not BCS S- 
wave pairing |jfl. The interaction leading to the for- 
mation of the singlet pairs appears not to be the tra- 
ditional lattice phonon mechanism underlying the BCS 
theory but a collective electronic interaction. The pair- 
ing gap has nodes in the k x -k y plane suggestive of D- 
wave hybridization in the 2-particle wavefunctions, and 
SC in the cuprates seems to be closely related to the an- 
tiferromagnetic (AF) insulator properties of their normal 
states. Furthermore, the formation of Cooper pairs and 
their condensation in high-T c compounds may be dis- 
tinct, with pair formation at a higher temperature scale 
than the formation of the SC state. This is reminiscent 
of grand unified theories in particle physics, where a hi- 
erarchy of symmetry breakings is implied by a Lie group 
structure that is broken spontaneously (and explicitly) 
on different temperature scales. Such observations ar- 
gue for a theory that is based on symmetries describing 
pairing more sophisticated than BCS and that can unify 
phases such as AF order and SC in the same theory. 

Dynamical symmetries having many of these charac- 
teristics have been developed extensively in nuclear struc- 
ture physics P,p|- There it has proven fruitful to ask 
the following questions: what are the relevant low-lying 
collective degrees of freedom, what are the microscopic 
operators that create and annihilate these modes, and 
what is the commutator algebra that they obey? It has 
been shown that dynamical symmetries having the rep- 
resentation structure described in Ref. Q] are realized to 
remarkably high accuracy in the spectrum and the wave- 
functions of large-scale numerical calculations Q] . 

Recently, S. C. Zhang and others have applied related 
ideas to high-temperature SC [p|-p[. In order to unify AF 
and SC order parameters, Zhang Tq] assembled these into 
a 5-dimensional vector order parameter and constructed 
an 50(5) group that rotates AF order into SC order. 
In this paper we proceed differently by closing a min- 
imal fermion algebra containing D-wave pairs of singlet 
spin (D-pairs) and triplet spin (7r-pairs), charge, and spin 



operators. Nevertheless, we shall recover an SO (5) sub- 
group of a more general U (4) symmetry. 

Thus, the recent discussion of Zhang's 50(5) symme- 
try applies directly to our results, but we extend this 
discussion in three ways: (1) The 50(5) subgroup is em- 
bedded in a larger fermion algebra that constrains the 
50(5) subgroup. (2) The SU(4) highest symmetry has 
subgroups in addition to 50(5) that may be relevant 
for AF and SC phases, and that aid in interpreting the 
50(5) symmetry. (3) We implement an exact dynamical 
symmetry using Casimir invariants of group chains and a 
corresponding collective subspace of low dimensionality. 

Let us introduce bilinear fermion generators: 

P12 = J29( k ) c U c -ki P12 = ^9*{k)c^ ki c k] 

k k 
k 

Qij = c k+Q,i Ck >3 ~ ^""^ C fc.t Cfc ij — 

k k 

where c\ i creates a fermion of momentum k and spin 
projection i,j — 1 or 2 = f or J,, Q = (7r,7r,7r) is an 
AF ordering vector, O is the lattice degeneracy, and 
g(k) = sgn(cosfcj; — cosk y ) with g(k + Q) = —g{k) and 
\g(k)\ = 1 (see Refs. |^,0|). The operator set (1) is closed 
under commutation, generating a U(4) Lie algebra with 
3 dynamical symmetry chains: 

D 50(4) x U(l) D SU(2) S x U{1) 
17(4) D 5*7(4) D 50(5) D SU(2) S x 17(1) (2) 
D SU(2) P x SU(2) S D SU(2) S x 17(1) 

that end in the subgroup SU (2) s x U (1) representing spin 
and charge conservation. The physical interpretation is 
aided by rewriting the generators of the U (4) algebra as 

Q+ = Qn + Q22 = X!( c fc+QT Cfc T + C k+Ql°kl) 
k 

, 5i2 + 521 . 5i2 — 521 5n — 522 \ 

5 = , -i , 
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si (Ql2 + Ql2 . Ql2 — Q21 Qll — Q22\ ,„% 

Q= { 2 2 ' 2 J (3) 

& = (5 fan - <&)> §(«u + ^22), -| (<Zi2 + 3ai)) 
^=(7ft)t £>t =p t 2 £) = pi2 M = ±(n-fi) 

where Q+ generates charge density waves, S is the spin 
operator, Q is the staggered magnetization, and n\ jf are 
the triplet D-wave pairs (Ref. g), the operators D\D 
are associated with singlet D-wave pairs, n is the electron 
number operator, and M is the charge operator. As we 
justify below, SU(2) p is associated with superconductiv- 
ity, 5*0(4) with AF order, and S0(5) with a "spin-glass" 
phase. 

The group SU (4) has a Casimir operator 

C su{i) = -7? + D^D + S-S + Q-Q + M(M - 4) (4) 

and the irreps may be labeled by 3 quantum numbers, 
(<7i, o"2, 03). We choose a model collective subspace 

|S> = Inxr^n.nrf) = (4)^(^(4)^(^10) (5) 

which is associated with irreps of the form (01, 02, 03) = 
(§,0,0). The corresponding expectation value of the 
SU(4) Casimir is a constant, (C au u\) — + 4). 

The full space is of dimension 2 2f2 but the collective 
subspace is much smaller, scaling as ~ f2 4 : 

Dim (§,0,0) = i(f + 1)(0 + 2) 2 (f +3) (6) 

This corresponds to truncation by a factor 10 54 for £1 ~ 
100. For small lattices we may enumerate all states of 
the collective subspace, yielding a simple model where 
observables can be calculated analytically. 
The charge density wave operator 

Q+ = X^ c !+QT c fci + c i+Ql Ck -i) = X^H'Xj ( 7 ) 

k rj 

(where the number operator for electrons on lat- 

tice site r with spin j) generates the U(l) factor in 
f/(4) — > U(l) x SU(A) and commutes with all genera- 
tors. Thus (S\Q+\S) — and charge-density waves are 
excluded from the present collective subspace. 

An ST/ (4) model Hamiltonian can be constructed from 
a linear combination of Casimir operators for all sub- 
groups: H = A + J2g ^GiCoi, where A and Ac { are 
parameters and the Casimir operators Cd are 

C so ( 5 ) =7T t -7r + 5-5, C so ( 4 ) = Q ■ Q + S ■ S 
C su(2)p = D^D + M(M - 1), C su(2)s =S-S 
C u{1) = M and M 2 (8) 

For conserved electron number the terms in M and M 2 
in Eq. (||) are constant. Since (C stl (4)) is a constant, 



by using Eq. (^) we can eliminate the 7?^ -7? term and by 
renormalizing the interaction strengths the SU(4) Hamil- 
tonian can be written 

H = H -G[(l-p)D^D+pQ.Q]+K efi S-S (9) 

with (1 — p)G = Ggg, pG = x cS an d K c ff as the effective 
interaction strengths, and < p < 1. The term Hq is 
a quadratic function of n and may be parameterized as 
T^o = £cff n + v c ff n(n — l)/2, where e e s and v e ff are the 
effective single-electron energy and the average two-body 
interaction in zero-order approximation, respectively. 

In the full SU (4) symmetry, D-wave pairing, antifcr- 
romagnetism, and 7r collective modes enter on an equal 
footing. The system has formed local SU (4) pairs, which 
fixes the length of vectors [SU (4) Casimir] but not their 
direction. Physically, the system is paired with SU(A) 
symmetry but is neither superconducting nor antiferro- 
magnetic since the energy difference in those directions 
are small on a scale set by the temperature. Neither AF 
nor SC order parameters have finite expectation values 
but a sum of their squares [the SU(4) Casimir] does. 

At this "unification" level, there is no distinction 
among these degrees of freedom, just as in the Standard 
Electroweak Theory the electromagnetic and weak inter- 
actions are unified above the intermediate vector boson 
mass. The full symmetry should hold while the temper- 
ature of the system is sufficiently high that anisotropic 
fluctuations in the directions associated with these col- 
lective degrees of freedom are negligible, but not so high 
that thermal fluctuations destroy the local ST/ (4) pairs. 
However, SU (4) symmetry is broken to its subgroups 
at lower temperature, leading to 3 dynamical symme- 
try chains with eigenstates labeled by the length of the 
SU(4) vectors, but with orientations no longer isotropic 
in the full SU(4) space. These low-temperature degrees 
of freedom may be interpreted as follows. 

The three dynamical symmetry limits SU(2), SO(4) 
and SO(5), correspond to the choices p = 0, 1, and 
1/2, respectively, in Eq. (||). The Hamiltonian, eigen- 
functions, energy spectrum and the corresponding quan- 
tum numbers of these symmetry limits are listed in Table 
I (where we introduce a doping parameter x that is re- 
lated to particle number and lattice degeneracy through 
x = 1 — n/Cl). The pairing gap A and the staggered 
magnetization Q, 

A={D^D)^ 2 g-(Q-Q) 1/2 , (10) 

may be used to characterize the states in these symmetry 
limits. As we shall see, each limit represents a different 
possible low-energy phase of the ST/ (4) system. 

(1) The dynamical symmetry chain ST/(4) D SU(2) p x 
SU(2) S D SU(2) S x U(l) corresponds to SC order and 
is the p — symmetry limit of Eq. (|^). The seniority 
quantum number v is the number of particles that do 
not form D-pairs. The ground state has v = 0, implying 
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that all electrons are singlet-paired. In addition, there 
exists a large pairing gap A = — x 2 ) 1 ! 2 , and Q = 0. 
Thus we propose that this state is a D-pair condensate, 
associated with the SC phase of the cuprates. 

(2) The chain SU(A) D 50(4) D SU(2) S corresponds 
to long-range AF order. This is the symmetry limit of 
Eq. (||) when p = 1. The 50(4) subgroup is locally 
isomorphic to SU(2)f x SU{2)g generated by 

F = ±{S + Q) G=\{S-Q), (11) 

where F and G are the total spin of electrons at even 
sites and odd sites, respectively. Therefore, the 50(4) 
Casimir operator can be expressed as 

O so(4) =2(F 2 + G 2 ). (12) 

The 50(4) representations can be labeled by the spin- 
like quantum numbers (F = w/2,G = w/2) and are of 
dimension (w + l) 2 . The ground state corresponds to 
lu = N and 5 = 0, and has n/2 spin-up electrons on 
the even sites (F = N/2) and n/2 electrons on odd sites 
with spin down (G = N/2), or vice versa. Thus it has 
maximal staggered magnetization 

Q = |fi(l -x) = \n (13) 

and a large energy gap due to the correlation energy Q-Q 
that inhibits electronic excitation, suggesting magnetic 
insulator properties at half filling. In addition, the pair- 
ing gap A = if2(x(l — x)) 1 / 2 is small (zero at half filling). 
Thus we propose that these states are associated with an 
AF phase of the cuprates. 

(3) The dynamical symmetry chain 5f7(4) D 50(5) D 
SU(2) S x U(l) corresponds to a phase with spin-glass 
character. This symmetry limit appears when p = 1/2 in 
Eq. (^) and has unusual behavior. At half filling, x = 0, 
the ground state is highly degenerate with respect to the 
number of 7r-pairs A, and mixing different numbers of tt- 
pairs costs no energy. The 7r-pairs must be responsible for 
the antifcrromagnetism in this phase, since only 7r-pairs 
carry spins. Thus the ground state in this symmetry limit 
has large-amplitude fluctuation in the AF order and we 
propose that this symmetry limit may be associated with 
a spin-glass phase of the cuprates. 

Away from the symmetry limits there are no exact so- 
lutions but the ground state properties can be studied 
easily using the coherent state method B. This is dis- 
cussed in a separate paper fiofl , but we quote one result 
of that study here to reinforce our point concerning the 
interpretation of the 50(5) symmetry. In Fig. 1, energy 
surfaces for the ground states for various particle num- 
ber n or doping x are plotted as a function of a quantity 
(3, which is related directly to the AF order parameter 
(see figure caption). Three plots are associated with the 
symmetry limits discussed above (p = 0,l,l/2) and one 



corresponds to a situation with a slight 50(5) symmetry 
breaking (p = 0.52). 

The energy minimum lies at (3 = for all the n's (or 
x's) for the SU(2) (Fig. la), and at j3 = [(1 - x)/4] 1 /2 
for the 50(4) limit (Fig. lb). In Fig. lc, there is a broad 
range of doping in which the 50(5) energy surface is 
almost flat in the parameter f3, implying large excursions 
in the AF order. Therefore, the 50(5) symmetry may 
be interpreted as a phase having much of the character 
of a spin glass for a range of particle numbers. Notice 
in Fig. Id that as doping varies the 50(5) Hamiltonian 
with slight symmetry breaking interpolates between AF 
order at half filling and SC order at smaller filling. Thus 
50(5) acts as a kind of doorway between SU(2) p and 
50(4) symmetries and this gives a precise meaning to 
the assertion || that the 50(5) symmetry rotates the 
AF and SC order parameters into each other. 

Our 5*7(4) model and Zhang's 50(5) model § have 
the same building blocks [the operator set (3)]. The es- 
sential difference is that we implement the full quantum 
dynamics (commutator algebra) of these operators ex- 
actly, while in Ref . |^] , the dynamics is implemented in an 
approximate manner: a subset of 10 of the operators acts 
as a rotation on the remaining 5 operators {D\D, Q], 
which are treated phenomenologically as 5 independent 
components of a vector. 

The embedding of 50(5) in our larger algebra has 
various physical consequences that do not appear if the 
50(5) subgroup is considered in isolation. We list three: 
(1) A phase transition from AF to SC at zero temperature 
and controlled by the doping emerges naturally, whereas 
in the 50(5) model a symmetry-breaking term propor- 
tional to a chemical potential has to be introduced by 
hand. (2) The present results suggest that the 50(5) 
subgroup is the appropriate description of the under- 
doped regime, but that the AF phases at half filling and 
the optimally doped superconductors are more simply 
described by our 50(4) and SU(2) p symmetries, respec- 
tively. (3) As we shall discuss in a separate paper |lj] , the 
present SU (4) theory leads naturally to the appearance 
of a pseudogap in the underdoped region [jllf , which oc- 
curs above the SC transition temperature T c and merges 
with T c near the optimal doping point. 

The temperature dependence of the phase diagram re- 
quires thermodynamic calculations that are in progress 
and will be presented separately. However, we may 
use the preceding discussion to construct the qualita- 
tive phase diagram illustrated in Fig. 2. First, Hq in 
the Hamiltonian ([)]) may be regarded as the energy scale 
for the U(4) — » U(l) x 5C/(4) symmetry, representing a 
fcrmion system in which electrons are all paired but with 
no distinction among the D-pairs and 7r-pairs. We may 
expect this symmetry to hold while the thermal energy 
is less than Hq. When the system is cooled, the pairing 
and AF correlations [H — Hq, sec Eq. (^)] become impor- 
tant, 5*7(4) breaks to its subgroups, and different low- 
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temperature phases will appear depending on the compe- 
tition between D^D and QQ interactions as a function of 
doping concentrations. From the preceding discussion, at 
zero temperature we expect the 5*0(4) AF phase to dom- 
inate at half filling, the SU(2) P SC phase to be favored at 
larger doping, and the intermediate doping region to be 
described naturally by the SO (5) spin-glass phase that 
interpolates between SC and AF behavior with doping. 
Thus, 51/(4) symmetry implies the schematic phase dia- 
gram of Fig. 2, independent of detailed calculations. 

In summary, an SU(4) model of High-T c SC has been 
proposed that contains three phases: A SC phase identi- 
fied with the SU(2) p dynamical symmetry, an AF phase 
identified with the 5*0(4) dynamical symmetry, and an 
50(5) phase extremely soft against AF fluctuations over 
a substantial doping fraction that we term (loosely) a 
spin-glass phase. Realistic systems may mix these sub- 
symmetries while retaining an approximate SU (4) sym- 
metry. Zero-temperature phase transitions are shown to 
be driven by the competition between the D-wave pairing 
and the staggered magnetization QQ interactions, as con- 
trolled microscopically by the hole-doping concentration. 
As we shall discuss in detail separately, this model leads 
naturally to the appearance of pseudogaps in the under- 
doped regime. Thus, we propose that high T c behavior of 
the cuprates results from a t/(4) symmetry realized dy- 
namically, and because this symmetry is microscopic its 
physical interpretation is accessible to calculation. This 
provides a possible understanding of the cuprate phase 
diagram and substantial insight concerning recent 50(5) 
theories of D-wave superconductivity. 
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FIG. 1. Coherent state energy surfaces for Eq. (9). 
The staggered magnetization Q is related to /3 by 
(Q) = 2D.f3 (n/2D. - Pl) 1/2 , where (3 is the value of (3 at 
the stable point, so /3 measures the AF order. 



FIG. 2. Schematic phase diagram for U (4) symmetry based 
on Fig. Id. The Ho in Eq. (9) is expressed in terms of hole 
doping x with x = 1 — n/Q.; e' cS = e c ff^i an d v^ ff =v c fff2 2 /2. 
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Table 1. The Hamiltonian, eigenstates and spectra in three dynamical symmetry limits of the SU(4) model. E g . s . is the 
ground state energy, AE the excitation energy, N = n/2 the pair number, x = 1 — n/fi, and k S0 4 = n c ff + Xeff- 



SU(2) limit: \ip(SU2)) = TV, «,S,m.) 


SO(4) limit: |i/j(S04)> = |JV, to, S, m s ) 




SO(5) limit: |i/>(S05)) = |t,JV, S,ro s ) 


<C s „(2 )p ) = i(n-«)(f2-« + 2) 


(C so( 4)> = tu(iu + 2), w = N-n 




<C so(5 > = t(t + 3), Q/2-t = JV-A 


H = H + k s S-S 


H = H + k so4 S-S 




H = H + k s S-S 


-G™ [c su(2)p - M(M - 1)] 






-G^ [C mW + M — C so(5) ] 


E t ... = H -\G^tf{l-x 2 ) 


B s .a. =H -i Xeff n 2 (l-x) 2 




E g .s. =H - ix cB V 2 (l-x) 2 


AE = i/G ( ^n + K eff S(S + 1), v = v/2 


AB = MX cff (l - z)^ + k so4 S(S + 1) 




AE = A iG^O + K cff S(S + 1) 


v = JV, AT - 1, . . .0; S = i/,i/- 2, ... or 1 


fi = JV, JV - 2, . . . or 1; S = iu, tu - 1, . 


. .0 


A = JV, JV — 1, . . . or 1; S = A, A - 2, ... 
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